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Abstract 

Time domain responses of porous media have been studied by some 
authors, but generally the possible descriptions have been given in the 
frequency domain. The aim of this paper, limited to materials with rigid 
skeleton considered as equivalent fluids, is to compare the descriptions 
by Johnson-Allard (J A) as well as by Pride-Lafarge (PL) with i) some 
analytical, approximate formulas, based upon asymptotic high frequency 
expansion ; ii) the exact formula by Zwikker and Kosten for the case of 
cylindrical pores. The paper starts with a short summary of the state- 
ment of the different general full frequency models (J A and PL). The 
Green function in the time domain is shown to exhibit interesting prop- 
erties of materials. In particular the maximum response depends on one 
dimensionless parameter only, which is denoted £ and is the ratio of the 
travelled distance to the product of the "frozen" sound speed and a charac- 
teristic viscous relaxation time. The distance £ is related to a time domain 
Stokes number. The numerical computation of the Green function is done 
by FFT, with some precautions, because of the importance of the higher 
frequencies on the response shape. The PL description is shown to be 
the best full frequency general model, but some discrepancies with the 
exact model appear at short times or short distances. When the distance 
£ increases from zero, the asymptotic expansion shows that the maximum 
of the Green function decreases first as l/£ 2 , then exponentially. 
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1 Introduction 



Long-wavelength wave propagation in homogeneous porous media is a classical 
topic. In particular for materials with a rigid skeleton considered as equivalent 
fluids, the classical homogenization theory gives us recipes to compute the com- 
plex, frequency-dependent, equivalent-fluid density and compressibility. This 
computation would require knowing, however, the detailed microgeometry of 
the solid. In absence of this knowledge, a widely used ID semi phenomenolog- 
ical frequency-domain model of propagation which depends on a small set of 
geometrical parameters of the structure is given by the well-known formula of 
Johnson et al. [1] for the density, and the Champoux-Allard [10] or Lafarge et 
al. [5] formulae for the compressibility. These approximate expressions (denoted 
J A) of the two constitutive functions essentially are the result of: 

1) a relatively precise description of the high-frequency limit, the two first 
leading terms of the two functions, density and compressibility, being exactly 
described in terms of the concepts of ideal-fluid tortuosity and characteristic 
lengths A (see Ref. [1]) and A' (see Ref. [10]), the viscous and thermal relaxation 
processes being "frozen" in this limit, 

2) a less precise description of the low-frequency limit, only the first non- 
trivial leading term of the two functions being exactly described in terms of the 
concepts of d.c. viscous and thermal permeabilities, fco (see Ref. [3]) and k' (see 
Refs. [9]|T2]), the viscous and thermal relaxation processes being "relaxed" in 
this limit, and finally, 

3) an assumption that the frozen and relaxed limits 1) and 2) are connected 
in the simplest authorized manner, i.e. (see pQ and [9,12]), by means of the 
simplest functions of frequency lo having their singularities and zeros lying on the 
imaginary axis (in the complex lo plane). Indeed, this assumption is equivalent 
to stating that the fluid velocity pattern at the pore scale is divergence-free, 
and this characteristics of the fluid motion tends to be verified (at least, as long 
as the geometry remains sufficiently simple) as a result of the long-wavelength 
separation condition which ensures that a homogenized medium exists. 

In this manner, a simple, resp. viscous and thermal, relaxation-transition 
description of the density and compressibility functions is obtained, that may be 
thought to be well- verified in a wide class of materials, as long as the wavelengths 
are large compared to the typical dimensions of the coarse-graining averaging 
volumes. 

Subsequently, the low-frequency relaxed limit 2) was made more precise 
by Pride et al. [TT| and Lafarge [12] (the next viscous and nontrivial thermal 
terms now being exactly described thanks to the introduction of the additional 
notions of d.c. viscous and thermal tortuosities ao and a' ), resulting in a 
slightly improved (denoted PL) description of the relaxation transitions of the 
two functions. 

In Ref. [13], Fellah et al. concentrated on the time-domain expression of the 
high-frequencies asymptotics implied by this model. 

Expanding, in the high-frequency limit, the PL density and compressibility 
in powers of the inverse Stokes number St^ 1 (defined as the ratio between 
boundary layer thickness and characteristic pore size), and retaining the (exact) 
zero and first order terms and the (model-dependent) second order terms, they 
derived an asymptotic time-domain pressure wave equation. Using fractional- 
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derivative and Laplace-transform calculus, they were able to rigorously solve 
this equation in instructive and elegant manner, through the calculation of the 
corresponding asymptotic Green function of the unbounded ID medium. Recall 
that a Green function or impulse response, as a function of time, extends and 
flattens when observed at fixed locations more and more remote from the spatial 
point where it originates. Like the wave equation used to construct it, however, 
the above asymptotic Green function was (partly) model dependent. 

Now, the results of the calculations made in [T3] seem to indicate that the 
terms of second order yield a significant effect on the Green function for resistive 
porous materials, and, in addition, that the effect is mainly an effect on the 
amplitude without noticeable distortion of the time wave pattern (see Figs. 2 
and 3 of paper (T3] ) . 

In the present paper, our purpose is threefold. 

First, paying attention to distinguish between the results that are general 
and the results that are model-dependent, we shall justify in another simple al- 
ternative analytic manner the finding of Ref. |13j summarized above. Indeed, by 
exploiting the model-independent fact that (with some restrictions on the geom- 
etry) the high-frequency asymptotic expansion of the wavenumber expands in 
successive powers of St -1 , and considering the o(2) expansion, we shall obtain a 
closed-form expression of a corresponding asymptotic time-domain Green func- 
tion G (2)(z,t). Exactly the same closed-form expression was obtained before 
by Polack et al. [TB] for the special case of cylindrical circular tubes. There- 
fore, like this special one, our general asymptotic o(2) Green function differs 
from the asymptotic o(l) one only by an amplitude factor exp(— C u x) indepen- 
dent of time: G (o\(x, t) — G m(x, t) exp(— G u x). This is the above-mentioned 
second-order attenuation with no distortion effect. 

Next, we hasten to point out that by construction, the models PL and J A 
unfortunately yield faulty values for the coefficient involved, and hence, 
faulty descriptions of this asymptotic second order attenuation effect. 

Finally, we try to assess in which circumstances this want of precision of 
the models PL and J A is manifested on the Green functions. To this end, we 
make detailed comparisons between the modeled and exact Green functions, 
computed by FFT, for the case of cylindrical circular pores. Recall that for 
this geometry, exact results exist since Zwikker and Kosten |14) . in frequency- 
domain, for the effective density and compressibility. It appears that the lack 
of precision of existing models is manifested mainly at sufficiently "short" prop- 
agation distances x (past the spatial origin x = of the impressed impulsive 
signal) and sufficiently "long" delayed times t—x/c (meaning the times t elapsed 
past the first arrival time x/c of the signal, with c = Cf / -Ja^ the frozen speed 
of sound in the medium and c/ the speed of sound in the fluid). In terms of the 
dimensionless distances £ = x/cQ and delayed times t = (t—x/c)/®, where 9 is 
a characteristic relaxation time (varying like the square of the pore size divided 
by the fluid kinematic viscosity), this means a limited region of the plane (£, 
t) near the origin. Outside this region i.e. at sufficiently "large" propagation 
distances x and/or delayed arrival times t — x/c (where "large" is inversely pro- 
portional, for a given shape of the pores, to the flow resistivity of the material), 
the existing models and especially PL, perform quite well. 

The paper is organized as follows. 

After the statement of the basic equations (section 2), the high-frequencies 
asymptotic wave equation and wave number are written (section 3 and 4), rc- 
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taining terms of order 0{2) (i.e. ~ St~ 2 ). As anticipated on general grounds, it 
is explicitly checked on the workable case of cylindrical circular pores, that the 
latter terms are not accurately described by the models J A and PL (section 
2). Using o(l) and o(2) asymptotic expansions of the wavenumber, two corre- 
sponding asymptotic Green functions in the time domain are derived (section 
4). These are the direct transposition of those given in [TS] for the case of cylin- 
drical circular tubes. By construction, our time domain o(2) Green function 
is obtained at the same order of accuracy as that derived in [T3]. The evalu- 
ation, by FFT, of all different time-domain Green's functions based on exact 
or approximate frequency-domain expressions, is discussed in section 5. The 
behaviors of the different expressions are compared and contrasted in section 6, 
in particular for the entirely workable case of cylindrical circular pores. 

2 Basic equations 

We start by recalling the form of the macroscopic Equivalent-fluid equations in 
the frequency domain (see [IWlll9|): 



where by definition, v and p are the macroscopic velocity and pressure obtained 
by coarse-graining (averaging) the microscopic fluid velocity and pressure fields, 
— iu) is the time derivative, pf and K a are the saturating-fluid density and 
adiabatic bulk modulus, and a(u>) and /3(w) are the dynamic tortuosity and the 
dynamic compressibility. Notice that for simplicity, isotropy or ID propagation 
along one principal axis is assumed, so that a{u) is a scalar. 

As every macroscopic theory, the Equivalent-fluid theory assumes a long- 
wavelength regime, i.e. a wide separation between the macroscopic wavelength 
and the typical porous medium lengthscale involved in obtaining the coarse- 
grained macroscopic variables. Now, the special form of Eqs. (TTJ) shows that 
the theory assumes much more than this. Indeed, as only the frequencies and not 
the wavenumbers are apparent in Eqs. (JT|), it appears that the Equivalent- fluid 
theory is certainly a simplified macroscopic description limited to describe the 
situations in which only the temporal dispersion - and not the spatial dispersion 
- effects are important [21] . Allowing for the spatial dispersion effects will lie 
outside the framework and scope of this paper. 

Coherent with this physical simplification which is well-justified in most ge- 
ometries (the condition to neglect spatial dispersion is certainly the same as that 
which allows us stating that the fluid velocity pattern tends to be divergence- 
free at the pore scale) an application of the homogenization technique [§] at 
order zero yields two different microscopic boundary value formal problems to 
be solved for determining the two constitutive functions a(ui) and P{lo). 

The first one describes the (velocity) response of a viscous incompressible 
fluid subject to an applied spatially uniform, time harmonic, bulk force source 
term: 



p f a(uj)iujvi = Vip ; K a 1 /3(u})icjp = V.u , 



(1) 




VII + V 2 w + e, 



(2) 



V • w 



(3) 
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in the pore space (with II a periodic or stationary random field, in periodic 
or stationary random geometries), and satisfying on the pore surface (no slip 
condition) , 

w = 0. (4) 

In this problem, 77 is the fluid viscosity, e is a dimensionless unit vector in 
the direction of the applied bulk force, and w is the scaled response velocity 
field (dimension of length 2 ). It determines the dynamic permeability k(u)) and 
dynamic tortuosity a(uj) introduced in the landmark paper by Johnson et al. pQ, 
by the relations: 

fcH = , ^ = (w> • e, (5) 
— iujpta{uj) 

where () is the coarse-graining averaging operation in the pore space, and <f> is 
the porosity. 

The other describes the (excess temperature) response of a thermal fluid 
subject to an applied spatially uniform, time harmonic, pressure source term: 

^l£L e = vH + i, (6) 
v 

in the pore space, and satisfying on the pore surface (no temperature-jump 
condition) , 

e = 0. (7) 

Here, 1 is a dimensionless unit constant representing the applied pressure, Pr 
is the Prandtl number, and 9 is the scaled response excess temperature field 
(dimension of length 2 ). It determines the functions k'(td) and a'(uj), thermal 
counterparts of functions k(ui) and a(ui) introduced by Lafarge 9,12 , and then, 
the effective compressibility /3(lj), by the following relationships: 

fc » = ■ J 14 * 77 \ - W ; /%) = 7 - 2 ^i. (8) 

2.1 Frozen and relaxed limits 

In absence of a complete information on the microgeometry, the two above- 
mentioned microscopic formal boundary- value problems cannot entirely be worked 
out and their exact detailed solutions (from which a(w) and /3(w) can in principle 
be extracted by coarse-graining, see Eqs.JSJ and flSJ)) are missing. Nevertheless, 
tacit, hidden, but very important assumptions on the geometry have been made 
in obtaining these problems (namely, that it is possible to neglect spatial disper- 
sion) , and these assumptions have had the consequence that in the first problem 
the stirring force e was set constant (and accordingly the velocity field v was 
represented by a divergence- free field), and in the second problem the stirring 
pressure term 1 was set constant (gradient- free). We shall later indicate what 
exactly mean these two characteristics of the obtained problems: they mean that 
the functions a(ui) and a'(u) have purely imaginary singularities and thus are 
relatively smooth functions on the real axis (this is property 3 mentioned in In- 
troduction). For the time being, we observe that in the limit of high-frequencies 
and low-frequencies, which we may call the frozen and relaxed, respectively, 
limits of the viscous and thermal relaxation processes involved, some general 
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characteristics of the solutions and corresponding functions a(ui) and may 
be sketched as follows. 



2.1.1 High frequencies (frozen limit) 

In the asymptotic frozen limit of high frequencies, the viscous and thermal 
terms V 2 w and V 2 # become negligibly small compared to the other terms. 
The fluid motions, except for vanishingly small viscous and thermal boundary 
layers at the pore walls, become close to those of an inviscid nonconducting 
fluid (77 = 0, k — 0, k is the thermal conductivity). Accordingly, the quantities 
iL ^ Pt w and l —— ' Pf 9 everywhere tend (except, of course, at the pore walls) to 
the "frozen" fields E and I verifying in the pore space, (with tp a periodic or 
stationary random field, in periodic or stationary random geometries), 

E = -V^ + e ; 1=1, (9) 
V E = 0, (10) 

and verifying at the pore walls (n is the normal on the latter), 

En = 0. (11) 

Now assume, following Johnson et al. [T] and Allard [5], that the pore-surface 
interface appears locally plane in this asymptotic high-frequency frozen limit. 
This is an assumption that the viscous and thermal boundary layer thicknesses 
8 = / and 6' — A / 2> ' p eventually become small compared to a character- 

istic radius of curvature of the pore surface. Then, the functions a(u)) and a'(u>) 
expand in integral power series of these thicknesses, which allows us writing a 
prior{]], 



1 / \ V 2 /1 

2aoo / V \ ( 1 



a(u) = aoo + =^[—^ +0. . 

A \—luipf J \—lU! 



A' \-iujpf Pr 

The intervening geometrical parameters ctoo and a'^ (dimensionless) and A 
and A' (dimension of length) must be some pore averages constructed with the 
frozen fields E and /. Indeed, detailed classic calculations made by Johnson 
and Allard [HUH], to which we refer the reader, show that they may be written 
as follows (see also [T^j and, in the most detailed manner for A, [S]): 

— $ 'i-n- 1 - < 13) 



1 Notice that the case of fractal geometry which modifies the exponent 1/2 in the first 
correction terms - see [J — is excluded by the assumption that the pore walls appear locally 
flat at the scale of the boundary layer thickness; the presence of sharp edges which modifies 
the exponent 1 in the second correction terms - see [5] — is also excluded by this assumption. 
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and, 

2 Is, V 2dS 2 fg^dS J Sp dS 

A ! v WdV ' A' J Vf PdV J Vf dV ' [ > 

where S p denotes the pore walls and Vf denotes the connected pore volume. 

Parameter A is an effective pore radius for dynamically connected pore sizes 
which was introduced by Johnson et al. [2] for the problem of electrical conduc- 
tion in the bulk fluid, perturbed by a thin, different conducting layer at the pore 
walls. Parameter A' is a length characterizing a simpler effective pore radius 
- twice the fluid-volume to fluid-surface ratio - sometimes referred to as the 
Kozeny-Carman radius; Allard |10j . inspired by Johnson's identification of the 
tortuosity high- frequency frozen limit, identified it as the thermal counterpart 
of parameter A. Mention that an erroneous expression for A was sometimes 
employed [3] ; the reasoning inaccuracjH was clarified and corrected in [5] along 
a line tentatively sketched in [6] (Appendix D). 



2.1.2 Low frequencies (relaxed limit) 

In the opposite relaxed limit of low frequencies, the viscous and thermal terms 
V 2 w and V 2 # eventually become much greater than the inertial terms '^ p/ w 

and _^Jp_fLt0 an( j the boundary layers extend to the whole fluid. Accordingly, 
the fields w and 9 everywhere tend to the d.c. - or "relaxed" - velocity and 
excess temperature fields Wo and 6>o verifying, in the pore space (with IIo a pe- 
riodic or stationary random field, in periodic or stationary random geometries), 

V 2 w = VUo-e ; V 2 = -1, (15) 
V-w = 0. (16) 

and verifying, at the pore walls, 

w = ; O = O. 

In this limit the functions a(to) and a'(uj) expand in Laurent's series, 



a(oj)= , # +ao + 0(-ioj) ; a'(cj)= . # p +a ' Q + {-iu>). (17) 
—lojpfko — iu>pfPrk 

The two first intervening geometrical parameters, ko (Darcy's viscous permeabil- 
ity) and k' (its thermal counterpart [9,12 ) on one hand, both having dimension 
of length 2 , and ao (viscous tortuosity) and a' (its thermal counterpart) on the 
other hand, both being dimensionless, are pore averages constructed with the 
relaxed fields Wp and 8q. Indeed, simple calculations show that they may be 



2 The ignorance, in the bulk, i.e. outside the viscous boundary layer or the perturbed 
conducting layer, of a perturbation contribution due to a perturbed ideal-fluid or electrical 
bulk flow field orthogonal to the leading bulk flow field E, and having, contrary to the latter, 
nonzero normal components at the pore walls. 
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written ( [7|H2]). 



^(w ) ■ e ; k' = (j>{9 ) , 



(18) 



a = 




,2 ' "0 — ,„ , 2 




(19) 



Mention that, among other things it is possible to show that, whatever the 
geometry, 



the equalities being satisfied only for the case of aligned cylindrical pores. 
2.2 Full-frequency models 

Let us now come back to our hypothesis that spatial dispersion is absent. This 
hypothesis is an assumption that the geometry is sufficiently simple not to have 
structures such as Hclmholtz resonators manifesting the presence of widely dif- 
ferent pore sizes (the resonator neck and cavity dimensions). It is in absence 
of such structures and effects that the long-wavelength separation condition en- 
sures that the microscopic flow-field may be considered divergence-free for the 
purpose of determining the microscopic fluid velocity pattern (see Eq. 
and likewise, the excess pressure field may be considered gradient-free for the 
purpose of determining the microscopic excess temperature pattern. 

Now, Johnson et al. have shown that, because the velocity field may be 
viewed locally divergence-free, the singularities of functions a(u) and k(ui) - 
poles, zeros, and branch points - necessarily are purely imaginary (see pQ Ap- 
pendix A). This also was shown by Avellaneda and Torquato [6j: assuming 
periodicity, a formal solution of the problem (J5J3J) directly expresses in tempo- 
ral domain, in terms of a distribution of relaxation times. 

For the functions a'(uj) and k'(uS), it was similarly shown that, because 
the excess pressure field may be considered gradient-free, the singularities of 
functions a'{uj) and fc'(w) also are purely imaginary (see [5] Appendix C). 

The wanted functions must therefore have very simple and smooth behaviors 
on the real axis. For example, calculations |12] inspired by the DRT formalism 
of Avellaneda and Torquato show that in the Laplace domain (s = —iuj > 0) 
the permeability function k(s) is always a strictly decreasing positive function 
on the real axis s > (and the same holds true in the same manner for the 
thermal permeability function fc'(s)). 

Finally, as suggested by Johnson et al., the functions fc(w) and a(w), as well 
as the functions k'(uj) and a'{uj), will be seeked as the simplest ones, satisfying 
both the frozen and relaxed limits and the condition that they have purely 
imaginary singularities. 

To proceed, let us introduce a Stokes number constructed using Johnson's 
concept of dynamically connected pore size A, viz: 



«0 > "oo > 1, 

k <k' ; a > a' ; A < A', 



(20) 
(21) 




(22) 
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Johnson has proposed the following expression of a{uj), 




a(u) = a 00 | I - ■ T -~ry:i -!- — Si* \ . ('->:-!) 

where 

M = ^ (24, 

is a dimensionless shape factor associated to the geometry. This expression is the 
simplest analytical ansatz that yields the exact first two terms at high frequen- 
cies (Eq. <jT2j) ) and the exact first leading term at low frequencies (Ea.(ljTT)l). and 
automatically satisfies the condition on singularities. Factor of 8 is introduced 
for convenience in Eq. (pM|) . so that M = 1 for cylindrical circular pores. 

Similarly, to describe the function a'(uj) let us introduce a second Stokes 
number corresponding to thermal effects: 



St' = A'J- lUJPfPl . 



(25) 



The ratio St' / St — y/Pr A'/ A is of the order of Pr 1 ^ 2 , not far from unity. 
Proceeding as did Johnson, Lafarge [12 proposed to write, 



C / )Lf/2 

^ = 1+ M^V 1 + 1^' 2 < (26) 



where 



M' - (27) 



<t>k r - 

is the thermal counterpart of shape factor M (M' = 1 for cylindrical circular 
pores), giving in turn a definite model for the dynamic compressibility 

W=7-(7-l)M«). (28) 

As one key recognition of this model - namely, that there is one length A' 
playing for thermal effects the role of Johnson's A for viscous effects - was due 
to Allard [TU], we refer to the combined modeling of functions a and (3 as to 
Johnson- Allard's (J A). 

Subsequently, Pride et al. [llj , studying oscillating viscous flow in convergent- 
divergent channels, evidenced the fact that Johnson's simple formula (|23l) (as 
was already anticipated in pQ) may significantly underestimate the imaginary 
part of permeability at low frequencies. To remedy this, they proposed different 
modifications of the formula. In essence these modifications produce formulae 
capable to account for the exact value of parameter ao, which, in such channels, 
may be significantly increased as compared to Johnson's. Notwithstanding, fo- 
cusing on the imaginary part of the permeability, they did not use and mention 
the parameter ao- Its identification (see Eq. (|19p) shows that it is a measure 
of "disorder" of the "Poiseuille" flow in the given microgeometry, so that it is 
increased not only by the convergent-divergent mechanism considered by Pride 
et al., but also, e.g., by the effect of a mere irregular distribution of solid inclu- 
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sions. Whatever the cause of the enhanced "Poiseuille disorder", it will make 
more significant the use of the Pride modified formula for the function a. The 
same general considerations hold true also, mutatis mutandis, for thermal effects 
and the function a' . 

Now, among the different modifications proposed in the first is the sim- 
plest one that is capable to yield the exact first two terms at high and low 
frequencies, and simultaneously, to automatically satisfy the condition on sin- 
gularities, whatever the values of parameters <j), fco, cyq, otoo and A (incidentally, 
this last nice feature of the formula is at variance to what was stated in this 
paper). Lafarge |12j expressed it in terms of the parameter ao and the same 
description was then immediately applicable, mutatis mutandis, to the thermal 
effects. The corresponding Pride-Lafarge's {PL) model formulae ar^l: 



a(u) = aoo 1 + — — [l-p +Px l + St 2 )), (29) 




"'M = i+lTTcJS \ 1 -P'+P'\l 1 + 7^^), (30) 



M' 2 

M'st' 2 \~ r 1 * y * ify/ 

= 7 _( 7 _i)/q/( w ), (31) 



where M and M' are as before, and p and p' are the new shape factors given 

by 

p 1 2k ^ ■ v > 1 2fc o (32) 

It must be realized that, most probably, as long as the geometry is relatively 
simple, and as a side result of the strong constraints imposed by the special 
location of the singularities of the functions, these PL expressions constructed 
using the information of one more exact term than JA's at low frequencies, 
will also and more importantly, describe in a more accurate manner all of the 
viscous and thermal relaxation. They reduce to JA's if one replaces (|3"2")) by the 
simplifications p = p' = 1. 

Now recall that in the frozen limit, as we have assumed a smooth pore-surface 
interface, our functions expand in successive powers of St -1 or St'^ 1 and the 
first two leading terms (O(0) and 0(1)) in these expansions are exactly described 
in terms of Johnson's and Allard's concepts of tortuosity oioo and characteristic 
lengths A and A'. We may try, using the improved but still approximate PL's 
Eqs. (|29M32p . to estimate the third 0(2) leading terms: 

<*» " -( 1 + l + ^ + (^)> _ < 33 » 

«"> = ^h- 1 >( 1 + i + l?^ 2 + (sF ! ))"- (M) 



3 Notice that in Ref. |13| there was a mistake of a factor 2 in the term under the root, 
without influence on further equations. 
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The second equation leads tcQ : 



«") = i + (7-i>(i; + ^-^ + °(^)). (35) 

Nevertheless, we should not hope that the above expansions yield anything 
precise as regards these third - model dependent - 0(2) terms. This is because 
the shape factors p and p' in the PL model (|29H30p are not fixed through the 
"frozen" conditions, physical in this context, that Eqs. ([33]) and (|35|) should 
match the exact frozen limits, but instead, through the physically very different 
"relaxed" conditions that the relaxed limits Eqs. (IT71) should be verified. The 
example of cylindrical circular tubes may serve to illustrate this in quantitative 
manner. 



2.3 The case of cylindrical pores 

For a material with cylindrical circular pores of identical radius R (say for 
simplicity, all parallel and aligned along the direction of propagation), the two 
boundary value problems, Eqs. J2H]) and (J2E]), determining a(w) and /3(w) are 
easily entirely stated and solved. In effect, these are nothing but the problems 
considered by Zwikker and Kosten |14j in simplifying (on account of the wide 
separation between wavelength and tube radius) the governing equations of the 
full Kirchhoff's theory of sound propagation in a cylindrical circular tube [15] . 
Thus, we may say that the conventional Equivalent-fluid theory (neglecting 
spatial dispersion and expressed by Eqs. (HE])) is the direct generalization, to 
the case of arbitrary geometry, of the classical Zwikker and Kosten's theory. 
Now, Zwikker and Kosten's result is that a(io) and /3(w) express via Bcsscl 
functions as follows: 

' l-x(w) 5 0(w) = l + (7-lM«Pr), (36) 



where x( w ) is the following relaxation function 



2Ji ((^i? 2 ) 1/2 ) 
XH = — V ',..„„ ' , N , (37) 



1/2 



(going from relaxed value 1 at low frequencies to frozen value at high fre- 
quencies). Using the known small-arguments and large-arguments series and 
asymptotic expansions of Bessel functions (or Kelvin functions), and comparing 
the results with those of the exact low-frequency and high-frequency limits (Eqs. 
(jTTJ) and fT2"]) '), the following parameters values are obtained: 

k = k' = (j)R 2 /8 ; a = a' = 4/3 ;a 00 = l ; A = A' = R. (38) 

^From Eqs. [[25]) and (|3"21) . the circular-tube shape factors M, p and M', p' are 
then identified as: 

M = M' = 1 ; p = p' = 3/4. (39) 



4 Notice that in Ref. [13] the last term —i/St' 2 in Eq J35t was omitted. 
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In order to see the extent with which these values may be used to estimate the 
frozen 0(2) terms for the case at hand (aligned cylindrical circular pores), let 
us write for this case the known exact asymptotic high-frequency expansions: 

«(«) - l + ! + ^+0(^); (40) 
m = ^-l)(i-^ + ^)). (41) 

If p and p' were fixed by the "frozen" conditions that the PL expansions ([33)1 
and (|35[) were to match the above exact expansions, they would be p — p' = 5/8. 
They are, however, fixed by other "relaxed" conditions, giving p = p' = 3/4. 
The difference (1/8) is not very large, but it produces a significant deviation 
on the 0(2) terms because of the factors of 8 in Eqs. (1331) and (1331) . Thus 
for instance, with the values (|3§)) . the term St~ 2 in Eq. ([35)1 has a coefficient 
equal to 2 instead of 3 in Eq. (|4T)]) . while the term St'~ 2 in Eq. (|3"5j) has a 
coefficient equal to —2(7 — 1) instead of —1(7 — 1) in Eq. (|4"Tjl ^l. These faulty 
PL coefficients, while significantly better than JA's (resp. instead of 3 and 
—4(7 — 1) instead of —1(7 — 1)), will obviously not yield accurate estimates of 
the exact 0(2) frozen terms. Here is a constitutive limitation of both models 
J A and PL that should be kept in mind. Its consequences on the time-domain 
Green functions will be illustrated later on. 



3 Asymptotic o(2) wave equation 



The ID wave equation that follows from Eqs. (1) is the following Hclmholtz 
equation: 

dx z K a 

Using the PL high-frequencies asymptotic limits ([35)1 and (1351) we get, 

C- 1 >(W-^ + 5Bf) +0 «' (43) 

and this is the form of the exact asymptotic result, except of course for the PL 
faulty values of p and p' , if we intend to use Eqs. (I32"j) . 

In the time domain, the corresponding asymptotic wave equation is written 
as follows: 

d 2 P (x, t) _ A d* P {x,t) _ B r* d 2 P (x,t')/dt' 2 dt , c dp(x,t) = Q (44) 



dx 2 dt 2 J VT^F dt 



5 In Ref. [13] because the term -4/St' 2 in Eq J35j was omitted, the total coefficient of the 
term in St'- 2 in Eq {35j was +2. 
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Comparison between (|42M43p and (|4"4"| shows that the coefficients are given by: 





A'VPr 



A 




(46) 



4A 



( 



8(1 -p') 



) 



A 2 1 



(47) 



+ 



4 + 



A /2 p r 



Again, this is the exact asymptotic o(2) wave equation in time domain, except 
for a faulty value of the coefficient C, if we intend to use Eqs. (|3"2"|) . These 
equations are the equations (14-16) of Fellah et al. [13], who computed the 
Green function for an infinite medium described by this wave equation (|44l) , by 
using the Laplace transform method^. We shall see in sections 4 and 5 that an 
alternative closed-form Green function may be obtained in simple manner for 
this same o(2) asymptotic frozen limit. Before detailing this result, let us give 
the general definition of the Green function we employ. 

4 A simple definition of a Green function and 
its asymptotic o(l) and o(2) frozen forms 

4.1 Definition; scaled form 

Another general method of defining and calculating a Green function in an 
infinite medium is by means of the effective frequency-dependent wavenumber 
k in this medium. Let us define our Green function G{x, t) (or impulse response) 
as a propagated Dirac delta impulsive signal S(t) imposed at x = 0, or more 
precisely, as the inverse Fourier transform of the propagation transfer function 
G ul (x,u) — enp(ikx): 



Setting to zero the viscosity and thermal conduction coefficients, no fre- 
quency dispersion arises. The wavenumber writes k = w/c where c, defined by 
the first Eq. (|45[) . is the frozen speed of sound c = Cf / '*Ja^. The Green func- 




tion (|4"5|) coincides with a Dirac delta propagated at this velocity c: G(x,t) — 



Setting to nonzero values the viscosity and thermal conduction coefficients 
the medium wavenumber k writes, 



making apparent a complex function h(St) that describes frequency dispersion 
induced by the viscous and thermal relaxation processes. The Green function 

6 Nevertheless, coherent with Eq. H35I I. an additonal term has been included in the bracket 
of the coefficient C in Eq. H46II . 




(48) 




5{t-x/c). 




(49) 
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PS]) writes, 



— oo 



2tt 



a; . 



f?(a;,f) = / — exp -iu(t - x/c) - i-h(St)x . (50) 



c 



This Green function g x (t — x/c) — G(x,t) now extends and flattens when 
observed at positions x more and more remote from the origin x = 0. Ex- 
plicit asymptotic expressions of the Green's function defined in this manner will 
now be obtained by considering high-frequency asymptotic expansions of the 
wavenumber k (meaning high-frequency asymptotic expansions of the function 
h). 

It will also be convenient to express these functions in scaled form, as func- 
tions of a dimensionless delay time r and a dimensionless position variable £. 
To take but one example, before being Fourier-transformed, the Green function 
G u (x, u) — exp(ikx) where k is - say - given by the model PL, may be viewed 
as a function of: 

(i) one dimensionless frequency variable, e.g. fl — w0 where O is the char- 
acteristic viscous relaxation time given by, 

e = ^, (5i) 

V 

(hence St 2 = -in), 

(ii) one dimensionless position variable, e.g. 

e-s. (B) 

(hence £ -1 / 2 = Av/ pfc/(xrf) can be regarded as a time domain Stokes number, 
when replacing — iu) by c/x in Eq. (|22p). 

(Hi) a number of dimensionless parameters characteristic of the form of the 
porous space but not of its absolute dimensions (aoo, M, p, M', p' , and A'/A), 
and, 

(iv) two dimensionless parameters characteristic of the fluid (7 and Pr). 

Suppose that the parameters (in) and (iv) of both the medium and the fluid 
are held constant. Function G bJ (x,uj) = exp(ifcir) then reduces to a function of 
n which is parametrized by £. There follows that the shape of the corresponding 
time-domain function will depend, in scaled form, on £ only, provided the time 
is counted in a dimensionless manner, e.g. for the time elapsed after the first 
arrival of the signal, 

r=(t-~)/e. (53) 

But function G(x, t) has the dimension of the inverse of time. Therefore, 

QG(x,t) = F s (£,T), (54) 

with F s a scale-invariant function that depends on the form of the pore space 
but not on its absolute dimensions: 

f°° dtt r , i 

F s (Z,t)= / — exp -iClr - ifi£ft(V=if2) . (55) 

J-oc 27T L J 
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4.2 Asymptotic o(l)'s and o(2)'s wavenumbers 

4.2.1 General asymptotic result and PL asymptotic result 

Recall that, since we assume a smooth pore-surface interface, in the frozen limit 
the product a(uj) /3 '(uj) / 'a oo expands in successive integral powers of St^ 1 (see 
footnote 1), hence giving the form: 



k = 



1 



a 



y/—iuj (—iu>) 
Constant c and B u are exactly given by: 



•o(2) 



1 

c 

2ni 



PfQ-oo 

K a ' 
1 + (7"1) 



nri 

B u =2. -n\ with 

V PS A 

A 



A'VPT 



(56) 

(57) 
(58) 



whereas C w is a new unknown constant. We have been unable to find any 
rigorous theoretical statements that can be made about this coefficient. In the 
framework of the PL model it is estimated as (see Eq. (|43p): 



Clc — 

2n 2 = 



V 1 



PS A 2 
8(1 -V) 
M 



ri2 with 



2( 7 -1)A , ( 7 -l)A 2 / 8(l-i/) 



AVPr 



A' 2 Pr 



M' 



3-7 



(59) 
.(60) 



However, for the reasons explained before, this is a rather loose determination. 
The relationships between these coefficients and the coefficients B and C are: 



B^ = B^c 2 /2 ■ C u = (Cc 2 - Bl)/2. 



(61) 



In what follows k a n\ and fc„(2) designate the asymptotic wavenumbers ob- 
tained by limiting the expansion of a(u})ff(u})/a ao to the first and second order 
exact frozen terms in St -1 : 



Mi) 



c 



B u 



Co(2) 



c 



1 



B u 



(■ 



y/—i(J 

4.2.2 The case of cylindrical pores 



(-iu) 



(62) 
(63) 



For the case of cylindrical pores, using Eqs. (j3"8")l and (l3l)l) . the expressions ([FT)) 
and (1591) of and CL become: 



# V />/ V VPr 



r 1 »? 1 



7 — 1 7" — 1 



/ Pr 



2Pr 



(64) 



As expected, the last expression is not the exact one. This was given by 
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Keefe [T7]: 



n 2 = 1 



7-1 7(7-1) 



/p r 



2Pr 



(65) 



For standard conditions in air, the exact result for 712 in C u is 1.08, while 
the approximated one is 0.29, i.e. more than three times smaller. Notice that 
Johnson- Allard's model would give a negative estimate of —1.27 for it. 

4.3 Frozen limit: Asymptotic o(l)'s and o(2)'s Green func- 
tions 

Taking for k its asymptotic o(l) expression (|62[) and making the inverse Fourier 
transform, it is a known result (20j that the corresponding analytical form of 
G(x,t) is, 



G (i)(x,t) — for t < x/c, 
G ^i)(x,t) — Fi(x,t) for t > x/c, 



where 

Fi(x,t) 



Q(x) 



^{t-x/cf/ 2 



exp 



Q\x) 



t — x/c 



, with Q(x) 



B^x/c 



(66) 



(67) 



The dimension of the function Q(x) is \fi. 

Likewise, taking for k its asymptotic o(2) expression (|6"5|) . it is also a known 
result of musical acoustics [TH] that the corresponding analytical form of G(x, t) 
is 



G a (2){x,t) = 0iort<x/c, 

G a (2)(x,t) — Fi{x 1 t)exp{—C u x/c)ioYt>x/c. 



(68) 



Thus, it is demonstrated that, as regards the Green's functions, the ef- 
fect of the second order terms on the inverse Stokes number (the effect of 
the term C u /(— iu) in the bracket (|63[1 ) is a pure multiplication by the fac- 
tor exp(— Cux/c): G ( 2 )(a;,i) = G ^(x, t) exp(— C^x/c). An expression ([5U)) is 
given for the coefficient C w , but, as we have argued, it will not be an accurate 
expression if we put JA's values (= 1) or PL's values (|52"]l of parameters p and 
p'. In general, the values p and p' that would give through ((59)) the exact value 
of parameter C u are unknown. They come from two new o(2) frozen geometrical 
informations (let's say, the values p^ and p'^ of parameters p and p' in Eqs. 
([55)) and (|M)) , such that these expansions are exact) that are simply not present 
in the parameters introduced so far (0, ctoo, A, A', k( h k' , ao, a' ). 

As examples of the scaling analysis, the formulae giving the o(l) and o(2) 
Green functions may be rewritten in the following scaled form: 



eG o(1) (x,t) 

QG oW (x,t) 



for t < 0, 
F al (t,r) forr > 0, 



(69) 
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where, 



1 nrf 



exp 



r 3/2 

where n\ is given by Eq. (|58[) and, 



, with 2ni = 1 + (7 - 1) 



A 

A'v/Pr' 



(70) 



8G ( 2 )(x,t) = 0forr<0, 

6G o(2) (a;, t) = F s i (£, r) exp(-n 2 £) for r > 0, 



(71) 



where n 2 is given by Eq. (|60|) . 



4.4 Relaxed limit: Asymptotic diffusive Green function 

The (relaxed) purely diffusive Green function Gdi / / can be obtained by assim- 
ilating /3(uj) with 7 and a(u>) with its first leading term r]4>/ (—iuipfko) in Eq. 
([T7j). In the relaxed d.c. approximation where /3 = 7 and a = T]<p/(— iujpfk ), 
one finds that the wavenumber expands as, 



k 



diff 



c \M8 y 



1/2 



(72) 



(we have used the relation (l2^|) 'l. Formally, this is the form (|r?2"|) with, in the first 
term, c replaced by 00, and, in the second term, B^/c replaced by the factor 



before yiuj. Straightforward calculations using the solution (|66ti67|) and the 
following dimensionless time f (I53[) with c = 00 and a prime to avoid confusion) 



then show that, 



t/e, 



&G diff {x,t) = 0forr'<0 
QG d iff(x,t) 



where, 



F.,diff(Sy) for r' >0, 
I /2t" C 



■ exp 



M r' 



,/tt V M r' 3 / 2 

5 FFT computations of Green's functions 



(73) 

(74) 
(75) 



In addition to the previous asymptotic analytic Green's functions, full Green's 
functions G(x, t) computed numerically by inverse FFT (Fast Fourier Trans- 
form) of the transfer function exp(ikx), can always be obtained as soon as exact 
or approximate full-frequency expressions of the wavenumber k in Eq. (|49[) arc 
known. These FFT computations are to be carefully done, however, as it is 
detailed below. 

For the case of arbitrary porous media, approximate models such as J A or 
PL may be used for k(w), leading to the computation by FFT of approximate 
Green functions G(x, t). For the case of simple workable geometries (e.g., cylin- 
drical circular aligned pores of identical radius) exact expressions are available 
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(e.g., Zwikker and Kosten's), leading to the computation by FFT of exact Green 
functions G(x,t). 

In section 6 we shall study particularly, in the case at hand of cylindrical 
circular pores, the relationships between the different, asymptotic or not, ap- 
proximate or exact, Green's functions. For later use, it will be convenient to 
distinguish and denote respectively by: 

G the exact (Zwikker and Kosten) Green function; 

Gja the full J A Green function (p = p' = 1); 

G PL the full PL Green function {p = p' = 3/4); 

G m PL the modified full PL Green function (defined by PL's expressions 
with p^p' = 5/8); 

An important point is that the diffusion solution Gdiff (Eq. ([74]) ) is used 
in our FFT computations to improve the accuracy of calculations. Indeed, in- 
stead of calculating directly the inverse Fourier transform G{x, t) of the function 
exp(ikx), we first compute the inverse Fourier transform G'(x, t) of the function 
exp(ikx) — exp(ikdif fx) (the latter difference being smoothly varying at low fre- 
quencies and rapidly decreasing at high frequencies), and next use the relation 
G{x,t) = G'{x,t) + G d iff (x,t), with G d iff(x,t) as given by the Eqs. (I73H75|) 
above (see Ref. [22]). The validity of the FFT computation has been checked 
for first order asympotic expression, with an accuracy better than 1%. 

The FFT computations can be compared with the following expressions: 

G e (x) the (frozen) o(l) Green function; 

G ( 2 ) the (frozen) o(2) Green function (p = p' = 5/8); 

GjAo{2) the J A o(2) Green function (p = p' = 1); 

Gplo{2) the PL o(2) Green function (p = p' = 3/4). 

With p = p' = 5/8, p = p' = 1, and p = p' = 3/4, the above expres- 
sion (jTTj) yields the Green functions G (2)(x, t), GjAo(2)(x,t), and Gpl (2) (%> t), 
respectively. 

Finally, the FFT computations can be compared to the solution of the 
asymptotic o(2) wave equation (|44[) . for both cases p = p' = 5/8 and p = 
p' = 3/4 (we choose to compute the solutions using FFT): 

Gweo(2) the (frozen) o(2) Green function (p = p' = 5/8); 

Gwepl (2) the PL o(2) Green function (p = p' = 3/4). 

These solutions being of o(2) as well as the solutions using the asymptotic 
o(2) wavenumber are expected to be very close to the latter. 

6 Computed results 

For the computation, the values of the parameters have been chosen close to 
those of material Ml in Ref. [13 . The porosity is 4> = 0.82, the flow resistivity 
a = 1960007Vm~ 4 s, the permeability fco = f]/cr = 9.225 10 -11 m 2 , corresponding 
to a radius R = 3 10~ 5 m. The temperature is 20°C, and the Prandtl number 
Pr = 0.71. The characteristic viscous relaxation time is O = 6.10 _5 s. 

6.1 Time responses 

The chosen length is 5.3mm (the dimensionless length is £ = 0.3152). Fig. 
[T] shows the results for the dimensionless Green function, obtained by FFT. 
The Pride-Lafarge description is compared with the exact Zwikker and Kosten 
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-0.5 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 

t/0 

Figure 1: Green function QG(x,t) with respect to time t/0 (both are dimen- 
sionless). (Black), solid line: Zwikker and Kosten formula (Eqs. (|36| and (l37l) . 
(Blue), dotted line Pride Lafarge description (Eqs. (l29l) to (|32jl ). (Green), thin, 
solid line: diffusive limit (Eqs. (|73)l and (|75p). For very long times, the diffusive 
limit is reached. 

formula. As expected, the full PL description is very satisfactory for short times 
(high frequencies) and long times (low frequencies). For very long times, both 
descriptions reach the diffusive (analytical) limit, i.e. the Poiseuille behavior is 
reached. 

In order to precise these results, Fig. [5] shows a zoom of the previous figure, 
and other approximations have been added. The modified PL Green function, 
denoted G m pL, can be compared to the PL Green function. As expected, it 
is more accurate for short time (during the signal rise) than the non modified 
PL function, but it is less accurate for long times, because the choice of the 
parameters p and p' has been done from the frozen limit expressions instead of 
the relaxed limit expressions. Otherwise the J A model yields less accurate 
results than both PL models. Notice that G (o\, the (frozen) o(2) Green 
function {p = p' = 5/8), which is a very simple analytical expression, leads 
to interesting results, except at long times. This will be confirmed by the 
calculation presented hereafter. 

6.2 Maximum values of the time response 

The previous results are concerned with a fixed value of the parameter £ , i.e. 
a fixed value of the thickness of the material layer. In order to compare the 
different descriptions for several values of £, we chose to compare the maximum 
values of the time responses (there is a unique maximum) . It is noticeable that 
for the asymptotic (frozen) expressions at the first order, F s i, the maximum 
value and the corresponding time r max are given by an analytical expression, 
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-0.5 1 1 1 1 1 ' 

0.3 0.35 0.4 0.45 0.5 0.55 

t/0 

Figure 2: Green function QG(x,t) with respect to time t/Q (zoom of figurc[T}. 
(Black), solid line: Zwikker and Kosten formula; (blue), dotted line: PL formula; 
(red), dashed line PL modified formula; (blue), mixed line: J A formula; (yellow), 
solid, pale line: frozen o(2)'s Green function G ( 2 ). 

obtained from Eq. d70| : 

Max(F sl ) = -= 
T _ ^ 2 

'max ^ "-IS • 

For the asymptotic expression of order 2, the time r max remains the same, 
while the maximum needs to be corrected from order 1 by the factor exp(— ri2^), 
according to Eq. (|TTj) . Fig. [3] shows the result for the exact Green function G 
and the two PL models, Gpl and G m pL, as well as the diffusive function Gdiff- 
It shows the neperian logarithm of the maximum value of F s with respect to 
the dimensionless space variable £. The two PL models seem to be very good; 
however we prefer to give more insight by subtracting the result corresponding to 
the exact Green function, as shown in Fig. As expected, the PL description 
is very good for high distances £, while the modified PL description is better 
for small distances £. Otherwise both are better than the J A description. In 
order to make possible a comparison with the asymptotic solutions (presented 
hereafter in Fig. [5]), the best asymptotic solution is shown in Fig. 0] It is bad for 
large distances. The transition range values of the distance £ is approximately 
between 0.036 and 1.7, corresponding to a range for the time domain Stokes 
number 5.3 > £~ 1/2 > 0.75. This range is similar to that accepted for the 
frequency domain Stokes number defined by Eq. (l22j) (see e.g. Ref. [17]). 

Fig. [5] shows that the PL description is slightly better than the first order 
of the frozen asymptotic, but the second order of the frozen asymptotic is the 



1-31" - 
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-(3/2) 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

i, 

Figure 3: Neperian logarithm of the Green function maximum. (Black) solid, 
thick line: exact Green function G. (Blue) dotted line: PL Green function 
Gpl- (Red), dashed line: modified Green function G m pL- (Green) solid, thin 
line: diffusive Green function, Gdiff- 

best approximation. Otherwise, for short distances £, the second order of the 
asymptotic wave equation WE exhibits the expected convergence to the results 
of the solution based on wavenumber expansion at second order, for both the 
coefficients obtained from the Zwikker and Kosten solution (p = p' = 5/8) and 
those of the PL model (p = p' = 3/4). Nevertheless, for higher distances, the 
asymptotic WE solution seems to be slightly less accurate than the asymptotic 
wavenumber solution: it is not easy to have an interpretation for this result. 

7 Conclusion 

For very short dimensionless distances, the value of the time maximum is in- 
versely proportional to the square of the distance, but when the distance in- 
creases, the decrease becomes exponential. This behavior is correctly represented 
by the second order of the frozen asymptotic. Nevertheless, if the coefficient of 
the second order is known for some geometries, such as for cylindrical pores, it 
is not known in general. The PL description remains the more satisfactory for 
the general case, but it gives faulty description of the exponential decreasing, 
because the coefficients p and p' are determined from the relaxed limit. The 
different figures presenting the results for dimensionless quantities can help for 
the choice of certain approximations. 

The second order solution of the wave equation, as presented in Ref. |13] 
differs (by definition) by the third order with the solution based on wavenumber 
expansion at second order, the latter being simpler to use in practice. Looking 
at the calculation made in the frequency domain, the figure 3 of this paper 
exhibits a ratio between the PL description and the J A one (with p = p' = 1) 
that is almost independent of frequency: this is confirmed by the asymptotic 
calculation of the wavenumber (Eq. (|71|l ). 

A problem of major interest is the use of the present investigation for the 
inverse problem, i.e. the determination of the parameters of a given material. 
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-0.3 1 1 1 1 1 1 1 1 ' 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Figure 4: Error on the neperian logarithm of the Green function maximum. The 
different solution are compared to the exact one, the error being the difference 
between the corresponding value and the eaxct value. (Blue) dotted line: Gpl- 
(Red), dashed line: modified PL G m pL- (Blue) mixed line: Gja- (Yellow) solid 
line: frozen o(2) Green function G ( 2 ). 

In this connection we have illustrated in this paper by direct FFT a simple 
fact. This fact is that the description of the time domain Green functions may 
be significantly improved by employing the full model expressions (instead of 
the asymptotic expressions). It suggests therefore that there is still a strong 
potential of improvement of the inverse methods of characterization based on 
recording transmitted and reflected pulses on different thicknesses of a material, 
provided the full expressions are exploited in the analysis. 

Assuming that we are dealing here with sufficiently well-defined and repro- 
ducible materials, we may suggest moreover, that this potential may be difficult 
to exploit fully, within the framework of current models including the most pre- 
cise PL. Indeed, the present work has highlighted one fundamental limitation 
of these models. It has shown that, even if it would be possible to gain precise 
experimental information on the unknown geometrical coefficients determining 
the high-frequency "second order" attenuation effect, it would not be possible 
to properly exploit this information in the framework of the current models. 
It is thus possible that working in the framework of a new model free from 
the highlighted limitation would present significant advantages for the precise 
inversion. These questions would have to be considered in further work. 
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